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FOREWORD 

This document presents the results of work performed 
by the Fluid Mechanics and Thermal Environment sections of 
the Lockheed-Huntsville Research & Engineering Center. This 
work was carried out under Contract NAS8-28057, "Study of 
High Altitude Plume Impingement," for the Aero-Astrodynamics 
Laboratory of Marshall Space Flight Center, The NASA Con- 
tracting Officer’s Representatives for this contract were Dr. 

T . F. Greenwood and Mr. D. C, Seymour, 
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SUMMARY 

A procedure is given for determining the radiation intensities in the 
base region of the space shuttle vehicle due to solid particle radiation ema- 
nating from the solid rocket motors of the shuttle. Results of an analysis of 
the Titan III and simulated solid rocket motor radiation intensities are pre- 
sented, A description of the gas particle flow fields of the Titan III nozzle and 
plume and a space shuttle solid rocket motor nozzle and plume is presented. 

A discussion of the gaseous Titan III flow fields is presented which utilizes 
the results of flow fields generated by a gaseous and two-phase method-of- 
characteristics computer programs. The development of a two-phase computer 
flowfield analysis program is discussed. An outflow correction theory is de- 
veloped which will be used to modify existing convection heat transfer methods 
for better heat transfer predictions on bodies immersed in rocket exhaust 
plume s .. 
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Section 1 
INTRODUCTION 

Successful design of the space shuttle vehicle will depend heavily upon 
the thermal environment produced by the exhaust of liquid and solid rocket 
motors which will be used for propulsion, attitude control and separation. 

The resultant thermal environment due to the rocket exhaust may be caused 
by heating due to. radiation from the plumes and/or direct plume impingement. 
Reliable analytic models must be available which will predict the environments 
in the exhaust plumes which contribute to the thermal environment encountered 
by the shuttle. 

The existing techniques for prediction plume impingement heating and 
radiation heating for liquid rocket motors have shown to be adequate during 
the Saturn program. However, solid rocket motors complicate the problem 
of predicting the environment within the exhaust plume. Solid rocket motors 
have metal additives in the propellant to increase the energy content of the 
system. These metal additives produce condensed solid and liquid particles 
in the exhaust gases which radically change the temperatures and pressures 
which are encountered within the plumes. Therefore, the techniques which 
were used to predict the plumes of liquid rocket engines must be revised and 
extended to handle the problems associated with solid propellant motors. 

Lockheed-Huntsville has reinitiated work on an equilibrium chemistry 
coupled two-phase flowfield analysis program under this contract. When this 
computer program is operational it will be able to describe the exhaust flow 
field of solid motors as well as existing liquid analyses describe the flow fields 
pf liquid rocket motors. Section 2 of this report briefly describes this computer 
program and lists its present capabilities. Ultimately the output from this 
program will be used for determining the environments produced due to plume 
impingement, solid particle radiation and gaseous radiation. 
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Section 3 presents the state-of-the-art techniques for predicting the 
thermal environment encountered by vehicles due to radiation from the solid 
and liquid particles which are present within a solid rocket motor plume. 
Because numerous vehicles in the past have been destroyed or severly 
damaged by this type of heating it is important that the analytic techniques 
for predicting this type of heating be well understood. Section 3 also pre- 
sents analysis and comparisons of the results of solid particle radiation 
studies performed for the Titan III and a simulated space shuttle solid 
rocket motor. 


When undertaking a study to advance the techniques in a field it is of 
interest to know what the differences in the results are between the old and 
new techniques. Section 4 presents the results of a study which predicts the 
exit plane properties of the Titan III 120-inch solid motor utilizing both solid 
and liquid rocket nozzle analysis programs. 

Of further interest of work performed during this contract are methods 
that were developed for improving the heating prediction techniques of the 
Lockheed-Huntsville Plume Impingement Computer Program (Ref. 1-1) de- 
veloped under Contracts NAS8-25511 and NAS8-21463. Current theories for 
predicting heat transfer rates on vehicles subjected to rocket exhaust impinge- 
ment do not utilize an adequate theory for predicting corrections whicn account 
for the outflow -induced thinning of the boundary layer. Section 5 presents an 
outflow correction theory that will be used at a future date to modify the existing 
convective heat transfer methods in the Plume Impingement program. 
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Section 2 

SUPERSONIC TWO-PHASE COMPUTER PROGRAM DEVELOPMENT 

The present concept of the space shuttle vehicle has numerous solid 
rocket motors for main propulsion, separation and abort capabilities. Metal 
additives are incorporated into the solid propellants of these motors to in- 
crease the energy content of the system and also suppress combustion pres- 
sure instabilities. The presence of these metal additives results in condensed 
products iii the exhaust gases which can result in several disadvantages . Since 
the condensed products can do no expansion work, their presence can only 
detract from the effectiveness of the nozzle to do expansion work in converting 
thermal to kinetic energy. Also, the presence of liquid or solid particles in 
the exhaust will contribute significant portions of radiation and plume impinge- 
ment heating on any structures which are either immersed or in close proximity 
to the exhaust plume. It is therefore important to know the physical properties 
throughout the nozzle and exhaust plumes. 

An extensive literature survey was made to determine what computer 
programs were available for calculating two -phase flow fields for use in the 
space shuttle base heating and plume impingement heating areas. Two basic 
types of tworphase computer programs were found to exist; (1) uncoupled two- 
phase programs such as Aeronutronics (Ref. 2-1) and; (2) coupled two-phase 
programs such as Kliegel's (Ref. 2-2). Uncoupled two-phase programs trace 
solid or liquid particles through an already generated gaseous flow Reid and 
do not consider the effect the particles have on the gas due to drag and heat 
transfer. Coupled two -phase programs include the effects of the particles 
on the gas as well as the gas on the particles. However, the existing coupled 
programs are primarily nozzle performance programs which do not include 
the flow chemistry {equilibrium or frozen) in the analysis. 
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The object of this study was to continue work done under Contract 
NAS8-21094 (Ref. 2-3), to develop a two-phase computer program which 
will describe two-dimensional and axisymmetric gas -particle flow fields, 
including equilibrium and frozen flowfield chemistry and the effect of shock 
waves. The two-phase capability has been incorporated into Lockheed's 
Streamline Normal Method-of-Characteristics Computer Program (Ref. 2-4). 

The governing equations have been rederived and incorporated into the 
two-phase streamline normal program. The equations were derived in a form 
in which the effects of equilibrium chemistry are included. The full develop- 
ment of the equations which will be incorporated into a future' document is 
based on the following assumptions: 

• The particles are spherically shaped. 

• The total mass of the mixture is constant. 

e The total energy of the mixture is constant. 

• The particle internal temperature is uniform. 

• The gas and particles exchange thermal energy by 
convection only. 

• The gas obeys the perfect gas law and is frozen or in 
chemical equilibrium, 

• The forces acting on the control volume are the pressure 
of the gas and the drag of the particles. 

• The gas is inviscid except for the drag it exerts on the 
•particles . 

• There are no particle interactions. 

• The volume occupied by the particles is negligible. 

• There is no mass exchange between the phases. 

• A discrete number of particles, each of different size 
or chemical species, is chosen to represent the actual 
continuous particle distribution. 

• The particles are inert. 

When the final checkout of the two-phasa program is completed it. will 
be able to operate in either the two-phase or gaseous-only mode, and the 


2-2 


LOCKHEED - HUNTSVILLE RESEARCH L ENGINEERING CENTER 



LtMSC-HREC TR D306389 


program will contain the same capabilities as Lockheed's Variable O/F 
Method-of-Characteristics program described in Ref. 2-5. The two-phase 
operation of the program will be able to handle up to ten different condensed 
species with each specie capable of having an input table of specific heat as 
a function of temperature. The equilibrium or frozen gaseous thermochemistry 
data will be calculated using a modified version of the NASA -Lewis Chemical 
Equilibrium Combustion Computer Program (Ref. 2-6), These data will be 
input to the two-phase program via tape or card input. The output from the 
two-phase program will contain all the gaseous and condensed species data 
which are necessary to, per form, any subsequent radiation heating or plume 
impingement heating analyses. Upon final checkout of the program a document 
will be published which contains a complete description of the program as well 
as a detailed input guide. 


2-3 


LOCKHEED- HUNTSVILLE RESEARCH & ENGINEERING CENTER 



LMSC-HREC TR D306389 


Section 3 

PARTICLE RADIATION ANALYSIS FOR THE 120 -INCH TITAN 
AND SIMULATED SHUTTLE MOTORS 

Design of the thermal protection system (TPS) for the space shuttle 
will require a thorough knowledge of the heating that results from particle 
radiation of solid motor exhaust plumes. Excessive heating to the base 
region has destroyed many of the earlier missiles during liftoff. Because 
of the weight penalty imposed, excessive' insulation for thebase' region is 
undesirable. It is therefore extremely important to be able to predict the 
heating to the base region incurred from the particle radiation. Much of 
the radiation from solid propellant rocket plumes is due to micron- size 
alumina particles. The complexity of the problem of predicting particle 
radiation requires knowledge of the following quantities: (1) particle environ- 
ment (temperature and pressure distributions of the plume); (2) the size dis- 
tribution of the alumina particles within the plume; (3) the scattering and 
absorption coefficients of the particles; (4) the index of refraction of the 
alumina as a function of temperature and wavelength; and (5) the radiance 
of the alumina particle cloud. 

Numerous methods are available for determining the radiance of a 
particle* cloud from a solid propellant motor but, currently, the most accept- 
able method is the Chrysler modified version of the Philco-Ford aeronutronic 
program (Ref. 3-1). This is an uncoupled program in that it considers only 
particle radiation with no gas radiation calculations performed. The purpose 
of the current study is to develop a confidence level in the program for future 
space shuttle studies by using the program to examine the particle radiation 
effects of two selected solid propellant rocket motors. 

The problem of predicting the particle radiation of plumes from solid 
propellant rocket motors can be divided into two disciplines. These are: 
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(1) particle environment and trajectory calculations, and {2) particle radiation 
predictions . Section 3 , 1 discusses the particle environment and trajectory cal- 
culations for the 120-inch Titan motor and a 1,3 scale version of the Titan 
(which closely approximates the shuttle engine). Relative comparisons are 
made between the particle trajectories and environments of the two motors. 
The particle radiation from the two solid propellant rocket exhaust plumes 
are discussed and compared. 

\ 

3.1 PARTICLE ENVIRONMENT AND TRAJECTORY 

Before particle trajectory or radiation calculations caii be performed, 
a gas flow field must be generated for the nozzle geometry of interest. The 
gas flow fields for the 120-inch Titan and the approximate shuttle engines 
were determined using the Lockheed -Huntsville Method-of-Characteristics 
(MOC) Computer Program (Ref. 2 -5) for the supersonic portion of the nozzles 
Real gas effects were considered in the flow field and were determined by the 
Chemical Equilibrium Composition Program (CEC) Ref. 2-6). The Aero- 
nutronics program assumes the gas flow to be ideal and one-dimensional from 
the rocket chamber to the supersonic section of the rocket nozzle. Problems 
therefore arose in matching these two theories, A point along the centerline 
of the nozzle downstream of the nozzle throat, at which there were sufficient 
characteristics data available, was selected for the interchange of the two 
theories. The static pressure, static temperature, Mach number and ratio 
of specific heats (gamma) were determined at this point from the MOC flow 
field. These local flowfield properties were used with a one -dimensional 
theory to compute new total conditions in the rocket combustion chamber. 

With the new total conditions, which were slightly different from the original 
chamber conditions, and gamma known, a flow field, assuming one -dimension 
theory, was obtained which matched the MOC flowfield properties at the given 
point downstream of the nozzle throat. The Aeronutronics program arranges 
the flowfield, generated by the two theories, by retting up a rectangular mech 
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of gas properties at various axial locations from the chamber to the nozzle 
exit plane and throughout the exhaust plume. This is done in order to pro- 
vide the necessary environments for the particles in the particle trajectory 
calculation. 

When the gas flow field for the nozzle and plume are completely speci- 
fied, subroutine TRAJEC, of the Aeronutronics program begins the particle 
trajectory calculations. The particle trajectory equations are solved by the 
method of finite difference. As always, the difficulty in using the finite differ- 
ence. .method lies in determining, the. time step such that the stability .require- 
ments are met for each of the equations to be solved. The basic equations to 
be solved by the finite difference method are; (1) particle location (x and y 
coordinates of the particle); (2) particle velocity; and (3) particle temperature. 
Reference 3-1 lists these equations along with four equations for determining 
time steps. The selection of the minimum value of the four time steps calcu- 
lated should meet the stability requirements in each of the equations. While 
attempting to determine particle trajectories for the 120-inch Titan it was 
found that the minimum time step calculated, from the four equations, was 
so small that the computer run time became excessively long. This minimum 
time step was to assure stability in the particle temperature calculation. It 
was decided not to use this time step but instead use the minimum value of 
the three remaining computed time steps. This was the correct decision in 
that the run time was greatly reduced with little or no effect on the particle 
temperature calculations. The only apparent result which could possibly be 
attributed to an instability in the finite difference equations were very slight 
temperature oscillations in the combustion chamber. These oscillations soon 
damped and were completely unnoticed as the particles enter the throat of the - 
nozzle. 

Figures 3-1 and 3-2 show the particle trajectories for the 120-inch 
Titan and shuttle motors, respectively. The particle trajectories for the . 
two motors do not vary much. This is an expected result however, since 
the shuttle engine and gas flow field are a 1.3 scale version of the 120 -inch 
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Titan engine and flow field. The chamber conditions for both engines are 
essentially the same. 

Particle temperature distributions for the Titan and simulated shuttle 
engines are presented graphically in Figs. 3-3 through 3-6. Figures 3-3 and 
3-4 show particle temperatures for centerline trajectories while Figs. 3-5 
and 3-6 show particle temperatures for the limiting trajectories. Particle 
temperature oscillations occurring along the limiting particle trajectories 
at the nozzle exit plane are evident in Figs. 3-5 and 3-6, There is an in- 
crease in the gas temperature due to the plume boundary shock which is 
attached to the nozzle lip during sea level operation. This leads to a gas 
temperature that is higher at the plume boundary, immediately downstream 
from the nozzle exit plane, than any other point in the plume. The gas tem- 
perature along this plume boundary remains higher than the temperatures of 
the plume core until the effects of the boundary shock diminishes. Figure 3-7 
is a radial temperature plot for the 120 -inch Titan plume. Temperature is 
plotted versus the non-dimensional parameter, R/R^, , for several X/R^, loca- 
tions, where R is the radial distance from the plume measured from the 
chamber and R^, is the nozzle throat radius. The temperature in the plume 
is essentially constant in the radial direction for each X/R^, location until the 
flow encounters the boundary shock represented by the discontinuities in the 
temperature plots. The low thermal capacitance of the smaller diameter 
particles, which are located on the outermost trajectory closest to the plume 
boundary, are extremely sensitive to the increased gas temperature. As a 
result of this sensitivity, the smaller particles experience a sudden. increase 
in temperature. The larger particles, on the other hand, have relatively 
larger thermal capacitances leading to a considerable time lag for the particles 
to adjust to the increased gas temperature. This allows the smaller particles 
to attain higher temperatures than the larger particles. As the effects of the 
shock wave diminishes downstream of the nozzle exit plane, the smaller 
particles cool as rapidly as they were heated and once again the larger particles 
have the higher temperatures. 
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Fig. 3-3 - Particle Temperature Plot for Centerline Trajectories (120 -inch Titan engine) 
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Fig* 3-6 - Particle Temperature Plots for Limiting TrajectorieS:{simulated shuttle engine) 
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3.2 PARTICLE RADIATION 

The general radiation equation for determining plume radiation heating 
incident upon a specific target of interest is 


Q = f f f l Hs 

A $ e 


sin6 cos9 d6 d<^ dX 


(3.1) 


At present, the computer program only calculates the spectral intensity, 

1x^0, emitted by the plume. The triple integration over the two angles, 0 
and and the wave length, X, is never performed and thus thei radiation 
heat transfer rate Q is not determined. The angles 0 and $ are input param- 
eters to the program and serve only to orient the line of sight from the plume 
to the target (see Figs. 3-8 and 3-9). 

Particle plume radiation intensity calculations are performed using the 
Mie scattering theory. The Mie scattering theory considers that the radiation 
which traverses an inhomogeneous medium such as an optically thin particle 
cloud is attenuated both by the scattering of the radiation into other directions 
and by the absorption of radiation by the particles. A digital computer pro- 
gram using the Mie scattering theory is used to determine the backs catte ring 
coefficient, (} , and the absorption coefficients as functions of particle radius, 
particle temperature and radiation wave length for the alumina particles in 
the rocket plume. At a given wave length, the scattering cross sections of a 
homogeneous spherical particle depend primarily on two parameters: (1) the 
ratio of particle projected area to radiation wave length, and (2) the index of 
refraction for the alumina particles. These two parameters must be obtained 
for input to the Mie scattering program before the scattering, o , and absorp- 
tion, o , coefficients can be determined. The latter of the two parameters, 
the complex index of refraction for the alumina particle, is the most difficult 
to obtain. These values must be determined by experimental means. The 
output of the Mie program (in the form of o , a , and p as functions of R , 

3L S P 

Tp and X) has been stored as permanent data for use by the radiation pro- 
gram. The experimental determined index of refractions were determined 
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Fig. 3-8 - Target Location and Lines -of-Sight Schematic 
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Fig. 3-9 - Schematic Showing End: View of Fig. 3-8 
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for wave lengths ranging from 0.5 to 5.0 microns and thus limits the prediction 
of the absorption and scattering coefficients. Reference' 3-1 contains the re- 
fractive index for both the molten and solid alumina state. 

A variety of different size particles exist at different states in the 
particle cloud of a two-phase rocket plume. The number distribution of 
particles in the plume is a function of position because the particle trajec- 
tories depend on particle size. The absorption and scattering cross sections 
are also functions of particle size. Since the particle size is a function of 
position within the plume the absorption and scattering cross sections are 
also "i func'ti on of position, ” 

The particle trajectory program computes the limiting and centerline 
trajectories and the particle properties along these trajectories for the various 
particle sizes in the plume. The particle trajectory data are written on mag- 
netic tape for use in the radiation program. Before the data can be used to 
compute radiation intensity for a specific line of sight through the plume , the 
trajectory data must be transformed into a more usable form. By performing 
a linear interpolation axially along each particle trajectory and then between 
the limiting and centerline trajectories a rectangular mesh of particle prop- 
erties whose grid size is entirely dependent upon the input value of the print 
interval, is obtained. The program checks on the line of sight and mesh line 
intersections to determine particle property variance through the plume. The 
line of sight through the plume is divided into a number of uniform e.qually 
spaced slabs. The number of these slabs is an input parameter. At each inter- 
face between successive slabs, along the line of sight, the computer program 
calculates the number density, particle size distribution, average radiative 
cross sections and particle energy flux. 

Initially, the size distribution in the chamber is given by the skew- 
symmetric formula 

* O p ] = 
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The size distribution is altered, however, as the particles move downstream 

by the influence of particle size on diffusion rates. The particle density ratio 

{the ratio of the local particle density to the chamber density, R^) is deter- 

mined for each particle size from the mass continuity relation. The particle 

density ratio, Rj, corrects the size distribution function, ^[R l t as the 

particles move downstream. At each slab interface the maximum particle 

radius, R max, is determined by surveying the y-coordinates of the limiting 
P , ■ 

trajectories, R max is broken down into a number of equal intervals, R , 

, P t P 

.such that, R , for instance, is equal to R max/50, Using R in Eq. {3,2) and 
P P P 

multiplying by the initial number of particles per unit volume in the chamber 

{input parameter RNTC) and the particle density ratio R^, the total number 
of particles per unit volume, N fc , present at the particular location is deter- 
mined by the following integration 



The scattering and absorption coefficients previously discussed, are 
dependent not only on wave length and temperature but also particle radius. 
The coefficients used in the radiation program at each slab interface are 
average values. Calculation of the average scattering coefficient, o , is 
determined from the following equation. 



(3.4) 


where N{R* ) is equal to ^(R ) • N • R,. The particle radiative properties 
o & and B are averaged over the size distribution that exists at each point 
along the line of sight as in Eq, (3.4) for o , Once all the radiative properties 
have been obtained for each slab interface on the particular line of sight through 
the plume, the radiation intensity equations are solved to determine the net 
spectral radiation intensity from the plume to the target. 
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The initial size distribution in the chamber for the Titan is shown in 
Fig. 3-10. This curve is of the form of Eq, (3.2) where a « 2.0 and b * 1,0. 

The initial size distribution in the chamber of the shuttle was assumed the 
same as the size distribution in the 120-inch Titan. 

There were two basic target locations for both the 120-inch Titan and 
shuttle motors investigated in the particle radiation analysis. Figure 3-8 
schematically shows two target locations and the corresponding lines of sight. 
This schematic can be applied to either the 120 -inch Titan or shuttle motor 
since the two targets are in basically the same location for each of the two 
engines. The two target locations for the 120 -inch Titan motor are: (1) 

X = 129 inches, Y = 180 inches, and (2) X = 209 inches, Y a 180 inches. The 
two target positions for the shuttle are; (1) X = 167.5 inches, Y = 234 inches, 
and (2) X = 272 inches, Y * 234 inches. The angles 8 ^ and 8 g are 300 and 
270 degrees, respectively, and are the same for the 120 -inch Titan and 
shuttle targets. The two angles (0^ and 82 ) along with the angle ^ (see Fig. 
3-9) orient the line of sight from the plume to the target. A radial survey 
of the plume was performed for each of the lines of sight for the two cases. 
This was done to determine the effects of varying the angle ^ on the radiation 
intensity calculations. Figures 3-11 through 3-14 are radiation intensity 
versus wave length plots for two targets in question for the 120 -inch Titan 
and shuttle engines. In particular. Figs, 3-11 and 3-12 show intensity as a 
function of wave length at the target 1 location for the Titan and shuttle engines 
respectively. Figures 3-13 and 3-14 are the same basic intensity plots for the 
120-inch Titan and shuttle engines determined at the target 2 location. 

Upper and lower bounds are placed on the intensity emitted from a 
particle plume to a specific target. The backs cattering coefficient, {3, deter- 
mines these bounds. The maximum radiation intensity emitted by the plume' 
(I Q (max)) incident on a target considers the backscattering coefficient to be 
zero, i.e., the medium the target sees has maximum emittance and no re- 
flectance. The other extreme, or the minimum radiation intensity emitted 
by the plume (I^fmin)), considers maximum scattering ( 3 = 1 ) i.e,, the medium 
exhibits maximum reflectance . I Q (j 3 ) is the radiation intensity value using the 
backscattering coefficient, |3, determined by the Mie scattering theory. 
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Fig. 3-13 r Target 2 Radiation Intensities for Titan 
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Fig. 3-14 - Target 2 Radiation Intensity for Space Shuttle 
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For a very thin medium, or a plume containing only a few particles 
which is usually the case in the small S-IVB ullage and retro motors (Ref. 3-1), 
the upper and lower bounds are nearly equal, whereas in a more dense plume 
with many particles the limits become wider. The exhaust plumes of the 
120-inch Titan and shuttle motors are very dense with particles due to the 
high concentration of alumina in the solid propellant for the two engines. The 
effect of the dense particle plumes are exhibited in the radiation intensity plots 

i 

shown in Figs. 3-11 through 3-4J.. Notice the wide spread (approximately two 
orders of magnitude) between the I^tmax) and I^fmin) plots in each of the 
figures. This variation indicates the importance of considering the back- 
scattering effects in the more dense particle plumes. In each of the figures 
(Figs . 3-1 1 through 3-14) notice the decrease in the spectral intensity curves 
at the 2 and 3 micron wave lengths. This decrease is due to the radiation 
blockage occurring as a result of the mean particle radius being on the same 
order of magnitude as the wave length. 

3.3 PLANK'S QUANTUM ARGUMENTS 


It has been shown by the quantum arguments of Planck and verified 
experimentally that for a black body the spectral distributions of hemispherical 
emissive power and radiant intensity in a vacuum are given as a function of 
absolute temperature and wave length by 


e* (x) = 


7TiX b (A) 


2 7TCj 


(e 2 - 1 ) 


(3.5) 


where (A) and i^ (A) are the emissive power and intensity, respectively, 
b v b 

and A is the wave length, T the absolute temperature and Cj and are 
constants determined by Planck. To better understand the implications of 
Eq. (3.5), it has been plotted in Fig. 3-15. Here the hemispherical spectral 
intensity is given as a function of wave length for several different values of 
the absolute temperature. One characteristic of Fig. 15 that is quite evident 
is that the intensity emitted at all wave lengths increases as the temperature 
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increases. Another characteristic of Fig. 3-15 is that the peak of the 
spectral intensity shifts toward the smaller wave length as the tempera- 
ture increases. The black body is defined as a perfect absorber and is 
also a perfect emitter. Its spectral intensity and therefore its spectral 
emissive power are only functions of the temperature of the black body. 

In each of the figures, there is a sudden increase in the spectral intensity 
values in the 3.0 to 5.0 micron wave length range. The spectral intensity 
value I^fmax) at the 5.0 micron wave length is comparable to that intensity 
emitted by a black body at a temperature between 2000° F and 3000°F (see 
Fig. 3-15). 

3.4 RESULTS 

Because the spectral intensity plots for the 120-inch Titan and shuttle 
motors reveal only slight variations in the spectral intensity values for the 
two motors, it can be assumed that the radiation heat transfer to the targets 
of interest is basically the same in each case. While it is true that at present 
there is no narrow view spectrometer data available for direct comparison 
with the analytical results, there is sufficient calorimeter data available for 
the 120-inch Titan motor. These data are presented in Ref. 3-2. Some of 

the calorimeter locations for the full-scale Titan motor static firings are 

t- 

shown in Fig. 3-15. Figure 3-16 also shows the X,Y, Z coordinates as well 
as the direction angles for each calorimeter. The radiation heat data was 
obtained from two 120-inch diameter motor static firings. Table 3-1 shows 
the results of one of these static firings. The principle targets selected in 
this analysis were the C012 and C009 targets (see Fig. 3-16). Since the 
simulated shuttle motor in this analysis is a 1.3 -scale version of the 120 -inch 
Titan motor, the Titan target coordinates were scaled by 1.3 to establish the 
shuttle target coordinates. This was done in order to maintain basically the 
same relative target positions between the two motors. Since the analysis 
indicated the same spectral intensity values at the two target positions for 
each motor, the experimental radiation heat transfer rate determined from 
the static firings of the 120-inch Titan motor can be assumed the same. 
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Fig. 3-16 - Calorimeter Locations and Coordinate System 
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Table 3-1 

RADIATION HEAT TRANSFER RATE PREDICTED 
BY EACH CALORIMETER 


Sensor 

No. 

q min 
/ Btu \ 

q max 
/ Btu \ 

q ave 
/ Btu \ 

q ave/F 
/ Btu \ 

\ft^-sec/ 

U 2 -s==; 

\ft^-sec/ 

\ft^-sec J 


11.7 

12.5 

12.0 

22.5 


10.7 

11.6 

11.1 

40.1 


5.0 

5.4 

5.2 

50.9 

C004 

11.5 

12,1 

11.7 

23.3 

C005 

11.8 

12.7 

12.2 

26.6 

C006 

10.5 

11.3 

10.8 

26.4 

C007 

5.0 

6.0 

5.2 

21.4 

C008 

8.6 

10.4 

9.0 

25.8 

C009 

8.7 

9.5 

9.1 

28.5 

C010 

0.9 

1.9 

1.2 

48.8 

C011 

2.1 

2.3 

2.2 

13.3 

CO 12 

2,7 

3.6 

3.3 

32.9 

CO 13 

0.7 

0.8 

0.7 

11.4 

COM 

— 

— 

— 

— 

0015 

i 
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Section 4 

COMPARISONS OF FLOW FIELDS GENERATED USING A GASEOUS 
METHOD-OF-CHARACTERISTICS PROGRAM AND A COUPLED 
TWO-PHASE FLOWFIELD PROGRAM 

One of the objectives of this contract is to determine how well the present 
liquid rocket flowfield anal/ sis programs predict the plume environments of 
solid rocket motors (SRMs). The Titan III 120 -inch SRM was chosen for this 
analysis. Two flowfield computer programs were used to analyze this motor. 
These programs are Lockheed'svMethodiof-Characteristics program (Ref. 2-5) 
(liquid or gaseous) and Kliegel's Two-Phase solid rocket motor analysis program 
(Ref. 2-2). For the purpose of comparing the results of these two programs, 
the nozzle exit plane Mach number and static pressure distributions were 
chosen. 

Using Lockheed's Method-of-Characteristics program, three different 
analyses of the Titan 120 -inch SRM were performed. The first two analyses 
consisted of running the nozzle using equilibrium chemistry thermodynamic 
data generated by the NASA- Lewis CEC program (Ref. 2-6) and ideal gas prop- 
erties obtained at the chamber conditions of the equilibrium chemistry. The 
third run was made using ideal gas chamber conditions with the specific heat 
rates modified to take into account the presence of. solid particles. The method 
used for modifying the specific heat ratio used here was obtained assuming no 
velocity and temperature lags between the gas and the particles, a gas gamma 
of 1.15 and a particle to gas mass flow rate ratio of 0,372. This method is 
described by Hoffman in Ref. 3-3, 

A fourth run was made using Kliegel's coupled two-phase nozzle program 
(Ref. 3-2). Kliegel's program is an ideal gas program. The gaseous thermo- 
chemistry data used in this run was the same as the ideal gas MOC runs with 
a gas gamma of 1.15. The particle to gas loading was 0.372, 
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The results of the four runs are presented in Figs, 4-1 and 4-2. Figure 
4-1 presents the exit plane distributions of static pressure. The gaseous 
method-of-characteristics run which matches the Kliegel's two-phase analysis 
is the run which was made using the modified specific heat ratio. The worst 
comparison is the equilibrium gaseous MOC analysis. Figure 4-2 presents 
the exit plane Mach number distributions from the four runs, The Mach 
number associated with the specific heat ratio MOC run is a gas Mach number 
modified in a similar fashion to the specific heat. The equations used are also 
presented in Ref, 3-3. 

The results of this study indicate that the best comparisons between 
gaseous analysis and a two -phase analysis are made using the modified specific 
heat ratio. The worst comparison is between the equilibrium chemistry gaseous 
run and the two-phase run. In this case the differences are due to both the two- 
phase losses and the differences between equilibrium chemistry and ideal gas 
assumptions. It will be very desirable to make this same comparison between 
programs when Lockheed's coupled two-phase program described in Section 2 
becomes operational with equilibrium chemistry. 
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Radial Distance from Nozzle Centerline at Exit Plane (in.) 


Fig, 4 


1 - Exit Plane Static PressureDistributions for 120 -Inch 
Titan III Solid Rocket Motor 
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Fig. 4-2 - Exit Plane Mach Number Distributions for 120 -Inch Titan 
III Solid Rocket Motor 
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Section 5 

PLUME IMPINGEMENT OUTFLOW HEATING CORRECTION 

Most current theories used in predicting convective heat transfer rates 
on shuttle-type geometries subjected to exhaust plume impingement do not 
utilize an adequate outflow correction theory. When a plume impinges on a 
surface at some angle of inclination, there is an outflow- induced thinning of 
the boundary layer. If current methods for predicting convective heat transfer 
without outflow correction theory-are used, the' predicted heating rates could 
be considerably lower than the experimentally measured heating rates. In 
many instances when sufficient experimental data are available, it is very 
difficult to make experimental-theoretical comparisons. Without outflow 
analysis to determine the severity of the thinning of the boundary layer and 
the corresponding increase in convective heating rates, erroneous experimental- 
theoretical comparisons may result. Convective heating rates determined by 
turbulent theory may compare more favorably with experimental data than 
laminar theory. This may lead to an assumption that the experimental data are 
turbulent which may. in actuality be laminar data. When the outflow correction 
theory is applied the laminar theory should match the laminar data and the 
turbulent -theory should match the corresponding turbulent data. 

The purpose of this study is to develop an outflow correction theory that 
will be used at a future date to modify the existing convective heat transfer 
methods on the Lockheed- Huntsville Plume Impingement Computer Program 
(PLIMP) (Ref. 1-1). 

5.1 THREE-DIMENSIONAL SURFACE FIT FOR PRESSURE DISTRIBUTIONS 

The calculation of inviscid surface streamlines and heating rates are 
dependent on the surface pressure distribution. Regardless of whether the 
pressure distribution is determined by experimental or theoretical methods, 

f 
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it must possess a certain characteristic. The pressure distribution must be 
such that it can be described by functional means. The function or functions 
used to describe the pressure distribution must also meet certain require- 
ments in that when differentiated they must yield continuous first and second 
derivatives in both the axial and circumferential directions. If the pressure 
distribution could be described by a discrete analytical function, this require- 
ment could be easily met. Unfortunately, however, this is not the case. In 
general, most if not all, pressure distributions cannot be described by an 
analytical expression, therefore a numerical technique must be used. The 
numerical method which seems to best fulfill the requirements is the Method 
of the Bicubic Piecewise Polynominal Functions. The cubic piecewise poly- 
nominal functions of one variable will be discussed first since their results 
are used to generate the bicubic polynominal functions. The theoretical pres- 
sure distributions used in this analysis were obtained from the Lockheed- 
Huntsville PLIMP program. 


5.2 CUBIC PIECEWISE POLYNOMIAL FUNCTIONS 


For functions of one variable, such as the cubic piecewise polynomial 

function, defines a function P(^) which assumes values Pj = P(0)j a *- a gi ven 

point 0. and matches the derivatives (dP/d^). at each ^.. See Fig. 5-1. 

1 J J 



Fig, 5-1 - Cubic Piecewise Polynomial Function of One Variable 
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Consider the following cubic polynomial 


. PjfoJ^Cj+C 2 ^ + C 3 ^ 2 +C 4 ^ (5.1) 

In order to determine the coefficients Cj - for the ploynomial Pj(^ ) such 
that when differentiated it will yield continuous first derivatives, a set of four 
simultaneous equations must be solved. Three of the four equations are of 
the form of Eq. (5.1) while the fourth equation is the derivative equation and 
is of the form of Eq. (5.2) 



= C* + 2C-tf . + 3 C . d* 2 
2 3 r j 4 r j 


(5.2) 



(5.3) 


By knowing the function values P.($ ) at <f>. where j = 1, 2 and 3 and the 

( 3P\ J J 

where j = 1 (calculated by Eq. (5.3), a set of four equations 


can be solved simultaneously for the coefficients Cj -C^. These coefficients 
will be used to compute any pressure or pressure-theta derivative in region I 
(see Fig. 5-1), The procedure for computing the coefficients for region I is 
repeated to determine the set of coefficients in each region for regions II 
through IV. The derivative where j =2, can now be calculated using 


the coefficients for region I and Eq. (5.2). Equation (5.1) where j =2, 3, 4 and 


the derivatiye 


ivati V e (If). 


make up a new set of equations which can be solved 


simultaneously for the coefficients in Region II. The end result is a set of 
coefficients in each region which when used in Eqs. (5,4) and (5.5) completely 
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specify the pressure and pressure derivative at any point ($). 

F(^) = Cj + C 2 tf + C 3 ^ 2 + C 4 tf 3 (5.4) 

|f = C 2 + 2C 3 ^ +3C 4 * 2 (5.5) 


The same procedure as has been previously discussed is used to determine 

the pressure and pressure derivatives in the x direction at constant 0 values. 

The cubic piecewise functions have been used to determine the pressure deriv- 
dP dP 

atives, and which are continuous derivatives at each grid intersection 
on Fig. 23. By knowing the pressure and the pressure derivatives with respect 
to x and y the mixed derivative may be determined by using the following 
equation 


a z p , i sV‘l-g(l) ff <»+»»- If (H 

d x 6 y ~ 2 ‘ Ay Ay 

The bicubic piecewise polynomial functions which describe the pressure 
surface in each region (Rj ^ - R.^) can now t> e determined by knowing the pressure 

- jp gp g 2p 

(P) and three derivatives ) at each grid intersection shown in 

Fig. 5-2. 

5.3 BICUBIC PIECEWISE POLYNOMIAL FUNCTIONS 
Consider the following bicubic polynomial. 


3,3 


p ( x * </> ) - C 00 + c in x + c ni^ + *“ c r< x $ 


10 


01 


33 


3 3 


= E E Cy xV j 


i = 0 j=0 


(5.7) 
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Examining the grid generated in Section 5.2 (Fig. 5-2) a statement can 
be made that applies in general to all bicubic polynomials. In each R.j there 
is one and only one bicubic polynomial (Eq. 5-7} which takes on specific values 

Qp ftp ft ftp 

of P, — , rr * t — 7TT at the four corners of the grid sections. This means 
that if these four quantities are specified at each of the four corners, of the 
grid section of interest, the bicubic polynomial function obtained will be unique 
to that grid section. When Eq. (5.7) is differentiated it will yield continuous fir st 
and second derivatives with respect to x and 0 at any point within the grid 
section. 

In each grid section, P is a bicubic polynomial (Eq. 5.7). Partial 
differentation of Eq. (5.7) with respect to 0 yields. 


8P, , . 

8 0 ~ 


3 

£ 

isO 


3 

£ 

i = i 


j C . . x 


*0 


j- 1 


with respect to x yields 


8 P 
8 x 


(x> 0) 


3 3 

E £ 

i=l j =0 


iC.. x 1 " 1 0 j 


(5.8) 


(5.9) 


and with respect to x and y yields 


8 2 P 
8 x 8y 


(x, 0 ) = ^ 2 


i = l j = l 


ii C.. x 




(5.10) 
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The second pressure derivatives with respect to x and with respect to 
y are present in Eqs. (5.11) and (5.12). 


8 2 P 


8 ^ 


2 s £ 22 j<j ' l) C V xl< 


;j - 2 


i= 0 j =2 


( 5 . 11 ) 


r? E E 

8 x i = 2 j =0 


i (i - 1)C.. x 


i-2 J 


(5.12) 


Each of Eqs. (5.7) through (5.10) is evaluated at each of four corners of the 
grid section of interest. These equations make up the set of 16 simultaneous 
equations to be solved for the 16 coefficients describing the surface in the 
particular grid section. 


To achieve a smooth continuous surface fit, as has been previously- 
discussed, continuity- in differentiation is a requirement. This quality of the 
surface fit is obtained by derivative matching (Eqs. (5.8) through (5.10)) at 
the interface between two successive grid sections as shown on Fig. 5-2. For 
example, let us examine regions R^ j and Rj., ^ -^6* 5-2. Suppose the coeffi- 
cients to be used in Eq. (5.7) to describe the surface over Region Rj^ have 

been calculated. When these coefficients are used in Eqs. (5.8) through (5.10), 

9 P 3P 9^p 

the equations must yield values for the and the g ■ at the inter- 

sections between Rj ^ and Rj^ that are identical to those previously calculated 
by the methods in Section 5.1. Likewise, the coefficients calculated for region 
Rj 2 when used in Eqs. (5.8) through (5.10) must yield derivatives which are 
also identical to those previously calculated at the interface between Rjj and 
Rj 2 * With this derivative matching continuous first, second, and mixed deriv- 
atives are assured. The end results obtained from the bicubic piecewise poly- 
nomial functions are presented in Fig. 5** 3. The pressure surface is completely 
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specified by individual sets of coefficients which uniquely describe a bicubic 
polynomial function in each region over the entire pressure surface. At any 
point on the surface , the pressure and pressure derivatives can be determined 
by using Eqs. (5.7) through (5.12) and the corresponding set of coefficients. 

5.4 CROSSFLOW EFFECTS BY STREAMLINE DIVERGENCE THEORY 

DeJarnette (Ref. 5-1) has performed extensive studies in the area of 
streamline divergence theory and its applications to convective heat transfer 
on reentry vehicles. His methods for determining inviscid surface stream- 
lines on these reentry vehicles will be adapted to problem area of determining 
the streamline patterns for particular shuttle type geometries subjected to 
exhaust plume impingement. 

To trace an inviscid surface streamline, an orthogonal coordinate system 
along the streamline must be determined. The geometry of the body is usually 
represented in a cylindrical coordinate system as r =f(x, 0 ) with unit vectors 
e^, e r and in the x, r and $ directions, respectively. (See Fig. 5-4). Since 
Euler's inviscid momentum equation (Eq. 5-20) only applies along a streamline, 
a coordinate system along the streamline which can be readily transformed into 
the original cylindrical coordinate system must be developed. Figure 5-4 shows 
the unit vectors in each of the two coordinate systems. 

The outer unit vector normal to the surface was determined by Dejarnette 

to be 


e^ = - sinT e x + cos V (cos 6 ^ e^ - sin 6^ e^) 


(5.-13) 


The unit vector perpendicular to the curve of the body but not necessarily 
normal to the body surface is 


e NN ~ cos e r ‘ sin6 <£ e 0 


(5,1'i) 
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Fig. 5-4 _ Unit Vector System Used to Determine Streamline Coordinates 
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By performing the vector cross products (Eqs. (5.14) and (5,15)) the 
unit vectors e^, (a vector tangent to the body surface) and e ^ ^ (a vector 
tangent to the body surface and perpendicular to the e^,) can be determined. 


A A 

e T = e x x e NN 


(5.15) 


A 

11 = 


A 



x 


(5.16) 


e . a unit vector in the streamline direction can be determined by the following 
s 

equation. 


A A A . * a A 

e s cos0 e n + $ui6 e™ 
s it l 


(5,17) 


eg defines a unit vector perpendicular to e g and tangent to the body surface 


A 

e 


p ~ 


A 

e 


A 

x e 


N 


(6.18) 


Performing the vector products and making the necessary substitutions 
the streamline coordinates in terms of the cylindrical coordinates can be 
written as - 


e = cos0 cosT e + (sin0 sinfi , + cosO cos 5, sinr) e 
s x p p r 


+ (sind cos 6, - cos0 s in 6 , sinF) e , 

0 0 0 


(5.19) 


e fl = - sin0 cosF e + (cos0 sinfi, - sinr sin0 cos 6,) ^ 
p x <p 0 r 


+ (cos 0 cos 6, + sin0 sinfi , sinl") e . 

0 0 0 


(5.20) 
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5.5 METHOD FOR DETERMINING THE STREAMLINE DIRECTION (6) 


Euler's inviscid momentum equation for determining the angle 9 (Fig. 25) 
along the streamline is 


DV v DV _ - VP 
DT " v DS " p 


(5.21) 


But V - Ve , therefore Eq, (5.20) becomes 
s 


£v = v £v g +v 2 = 

DT V DS s tv DS p 


'(5.22) 


In using Euler's equations, it is assumed that there is a constant entropy level 
on the body surface. 


Equation (5.21), expressed in streamline coordinates, becomes 


D9 

DS 


sinT 


Da 

DS 


pv‘ 


± 

h 


3P 

ap 


(5.23) 


where h is the metric coefficient that will be discussed later in the report. 


DeJarnette, in his analysis, determine the transfomation operators re- 
lating the streamline coordinate system to the cylindrical coordinate system 
(Eqs. (5.24) and (5.25)). 


1 _a_ x ax j*. i M a. 

h 8(J ■ h 3(5 3x h df ~B<j> 
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A 

= e 

x 


JL 

P dX 


V e & 

f 


_a_ 


g 

= - sin9 cosr -5— 
8 x 


(cos0 cos 6 , + sin0 sinS.sinD 


8 

Bth 


(5.24) 


and 


D J 8_ 1 8x 1 . 1M1 

DS - h s ae ~ h s ae ax + h g ae a0 


A 

= e 

x 


% JL + ‘i ' *« J. 

s ax f a$ 


= cos0 cosF 


a 

dx 


(sin0 cos6, - cos0 sin 5, sinT) 


80 


(5.25) 


Using the transformation operators (Eqs. (5.24) and (5.25))» Eq. (5.23) 
can be expressed in cylindrical coordinates in the following manner. 



- sin0 cosT 



+ 


(cos0 cos 6, + sin0 sing, sinF) 
0 o 

i 



- sinr 


cos0 cosF 


ap 1 

8x 


(sin0 cos { 


- cos0 sin6, sinT) Q 
0 _9a 

f 80 


(5,26) 
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Also using Eqs, (5.24) and (5.25), the geometric location derivatives 
with respect to streamline direction were determined to be 


D<j 

DS 


Dx 

DS 


=' cos© cosr 


sin0 cosfi . - cos6 sin5, sinT 


(5.27) 


’(5:28) 


- 3r- r A i ■ ’ 


For a cylinder, which is the geometry under consideration in this anslysis, 
the angles Tand 6^ are identically equal to zero reducing Eqs. (5.26) through 
(5.28) to 


* (72) [- sine ■& (p/p.) * S 2 r- if (p/p,)] 


(5.29) 


Dx 

DS 


= COS 9 


(5.30) 


D6 sine 
DS " f 


(5.31) 


Equation (5.29) can be solved by numerical integration, once the pressure deriv- 


atives 


a<p/pj 3 (p/p s ) 

and , - " are determined. These derivatives are deter- 


dx “““ dtf, 

mined by differentiating the piecewise bicubic polynomial functions describing 
the surface pressure distribution on the cylinder at the correct geometric 
location determined by Eqs, (5.39) and (5.31). Once the initial streamline direc- 
tion (6) is assumed, Eqs. (5.30) and (5.3>1) can be numerically integrated to deter- 
mine the correct geometric location (X, 0 ) at which the bicubic functions are 

& (p/p ) a (p/p ) 


evaluated for the 


3x 


and the 


a 0 


, in Eq. (5. 29). 
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5.6 METHOD OF DETERMINING THE METRIC COEFFICIENT, h 

Before presenting the methods for computing the metric coefficient, h, 
perhaps the physical meaning of the metric coefficient should be presented. 



Fig. 5-5 - General Orthogonal Curvilinear Coordinate System 

For any general orthogonal curvilinear coordinate, say £, (3 on the surface 
of a body, the distance between adjacent lines at constant £ or (3 will vary 
from point to point. If increments in these quantities are to be related to 
actual, physical distances along the surface (dS and df), it is therefore 
necessary to introduce metric coefficients say h. and h > such that h, df} 
and h^ d£ are, respectively, small distances (dS and df) measured in the 
| and (3 directions (see Fig. 5-5). Since the divergence between successive 
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streamlines is determine by h, hj is not of particular importance. The term 
h dP is equal to the integration step size dS for solution of the equations for 
h. 

t 

1 = 

Dejarnette (Ref. 5-1) suggests^the. following. differential equation for ihe 
metric coefficient, h, 


Dh 

DS 



sinr 



(5,32) 


Using the various substitutions presented in Ref. 5-1, Eq, (5-31) can be 

30 

expressed as follows to eliminate ypr, 



s 

P V 2 
*s «> 


2 6 

S2. i. d (p/p ) ] (3 - 

P yZ h 9 p tP/ VJ 


M 2 ) 


+ 


P 

s 


P V 2 


£s 

p 



1 JL 

h 9 p 



+ 


2 

cos 


r cos 6 
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4 


9 F 9cr _ 9a 9 P 
9x *90 9x 90 


(5.33.) 


Equation (3.32) is a generalized equation which applies to any geometry. 
Reducing this equation for a cylinder yields 
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X P 2h 

h DS 2 


s / . a 9P , cos 6 8P \ 

pV *A dZ + — 8?j 


(3 - nT> 


, 1 Dh,7 P 3\f. nnfl 9P sine 9 P\ 
- <h^>^2/\ COfle ©^ w) 

+ P S /l . 

+ 2 \T 31 

pV V 


in© cos© ■§? + sin 2 © — — ^ 
9< * 9 x 2 


' - -2/f sin© cos© + “ 


2 ‘ cos 2 © 9 2 P^ 


f 2 ‘ d<f> 


(5 


9 P 9F 

The derivatives, -r — , -r~r » 
dx d<p 


9 2 P 

9x 2 


2 2 
9 P 9 P 

— — and the ■ , are determined from 

9* 2 9x8* 


the surface pressure fit. This equation along with Eqs. (5.29) through (5.31) 
are solved simultaneously by a numerical integration routine (Runge-Kutta) 
for the geometric location of the streamline (x,^), the streamline angle (0) 
and the metric . coefficient, H. 


5.7 RELATIONSHIP OF THE METRIC COEFFICIENTS TO CONVECTIVE 
HEAT TRANSFER THEORY 

5.7.1 Laminar Convective Heat Transfer 


Heat transfer through a laminar boundary layer is calculated using the 
integral form of the energy equation. The effects of variable freestream 
properties (velocity, density and pressure) are accounted for through use of 
the appropriate transforms of the flat plate solution. Non-constant properties 
through the boundary layer are also accounted for. 

The effects of variable gas properties through the boundary layer on 
laminar boundary layer heat transfer are determined by evaluating the gas 
properties at a "reference enthalpy" and using these values in the constant 
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property solutions obtained by Blasius, Based on this method, the convective 
heat transfer rate to the wall is evaluated using 


q = g, (H - h ) 
^ s h ' r w 


(5.35) 


where the adiabatic wall or recovery enthalpy is 




e „2 


H = 


C dT + -=^r- 
P 2 Jg 


W'V': r ’ ^b. ^ ’ 

(5.36) 


where 


= Vpt" 


(5.37) 


The enthalpy conductance is evaluated from the Stanton number (St) 


St = 


s h Nu 


p* V 


* * 
e Re Pr 


(5.38) 


The "reference enthalpy" h can be determined by 


h* = h + 0. 5 (H - h ) + 0.22 (H - hj 

6 S 6 r c 


(5.39) 


The Blasius flat plate solution assuming constant property flow is 


Nu 

VRe 


= 0.332 Pr 


1/3 


(5.40) 
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By making the necessary subsitutions and noting that 


* P* V. X, 

Re s — 


e 


-Ji* 


(5.41) 


Equation (5,34) becomes 


• _ 0.332 
q = Pr* 2/S 


(p V V e ) 

3”7T 


0.5 


(H - h ) 
r w 


(5.42) 


The term Xj^ in Eq. (5.42) is the characteristic running length and can 
be obtained by numerically integrating the following equation along a flowfield 
streamline 


S 

X L = » 1 2 f h 2 dS (5.43) 

^ p liVh* J 6 

e o 

where S is the distance along the streamline and h is the metric coefficient. 
5.^.2 Turbulent Convective Heat Transfer 

The equation for turbulent convection heat transfer is 


q = 


0.0236 (p * y , 0.8 

6 


Pr 


* 0.66 


© 


<H r - h w ) 


(5.44) 


where X^. is the turbulent characteristic running length and can be calculated 
by Eq. (5.45) ' 


X, f = 


* * 

P M V, 


1,25 


/ 


* * 12 ^ 

p a v h i,Z5 ds 


(5.45) 
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Again, as in the laminar case, the metric coefficients contributions to the 
characteristic running length are revealed. 

5,8 RESULTS 

A digital computer program using the theory presented in the preceding 
section was developed to calculate in viscid surface streamlines on a cylinder 
subjected to exhaust plume impingement. The pressure distribution used by 
the program was calculated by the Lockheed -Huntsville PLIMP program 
* -(Ref ; 2-3) and is presented graphically in KFigv 5-6 . • Figure- 5-7 is a plot of- 
the inviscid streamline coordinates ( <j> versus x) for several initial streamline 
starting angles, 0, while Fig. 5-8 is a plot of the metric coefficient versus 
distance X for one of these streamlines (0 - 0,10). Notice in Fig. 5-7 how 
the streamlines diverge at an increasing rate as they proceed axially down 
the cylinder. This divergence, i.e., the spreading of the streamlines, is 
reflected by the metric coefficient plot in Fig. 5-8. The metric coefficient 
increases at a moderate rate to X = 4 inches. From X = 4 inches to X =4.5 
inches, there is a sudden abrupt increase in the metric coefficient that can 
be attributed to the increased spreading between streamlines. 
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Fig. S*6 - Surface Pressure Distribution for a cjylind 


0,32 


LMSC-HREC TR D3Q6389 



5-22 





LMSC-HREC TR D 306389 



0 1 2 3 4 5 6 

Axial Length Down Cylinder, X (in.) 


Fig. 5-8 - Metric Coefficient vs Axial Length Down Cylinder 
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Section 6 
: CONCLUSIONS 

Development of an equilibrium coupled two-phase computer program 
for solid rocket motor nozzle and plume analysis has been initiated and 
described. This program will provide a description of the exhaust of solid 
rocket motors which may be utilized in performing base heating and plume 
impingement analysis for the space shuttle vehicle. A comparison of tech- 
niques for predicting solid motor exhaust flowfields was presented which 
indicates significant differences in the predicted flow fields obtained using 
a two -phase computer program and a gaseous flowfield program. 

A plume impingement convective heat transfer theory with outflow 
corrections was presented that could be incorporated into the heat transfer 
portion of the Plume Impingement program which would result in an outflow 
corrected heat transfer analysis which may be important for shuttle-type 
vehicle applications. 

The solid particle radiation analysis for the Titan 120 -inch motor and 
simulated shuttle motor presented in this report is based on several quantities 
including: the particle environment, particle size distribution, the scattering 
and absorption coefficients of the particles, the index of refraction of the 
particles as a function of temperature and wave length and the radiance of 
the particle cloud. The technique used in this analysis has shown that better 
data on scattering and absorption coefficients and index of refraction of the 
alumina particles is desirable for any future analysis of shuttle motors. 

Also a technique for more easily obtaining radiation data once more than one 
line of sight should be considered. With the development of the equilibrium 
coupled two- phase analysis a more reliable flow field for obtaining radiation 
heat transfer data will be realized. 
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